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Abstract 

We have studied the three quark systems in an Effective Hamiltonian approach in QCD. 
With only two parameters: the string tension a = 0.15 GeV 2 and the strong coupling 
constant a s = 0.39 we obtain a good description of the ground state light and heavy 
baryons. The prediction of masses of the doubly heavy baryons not discovered yet are 
also given. In particular, a mass of 3637 MeV for the lightest ecu baryon is found by 
employing the hyperspherical formalism to the three quark confining potential with the 
string junction. 



The discovery of the B c meson |IJ demonstrates that new sectors of hadron physics are 
becoming accessible to experiment. In particular, the existence of doubly heavy baryons 
is a natural consequence of the quark model, and it would be surprising if they did not 
exist. Data from the BaBar and Belle collaborations at the SLAC and KEK B-factories 
would be good places to look for doubly charmed baryons. Recently the SELEX, the 
charm hadroproduction experiment at Fermilab, reported a narrow state at 3519 ±1 MeV 
decaying in A^K~n + , consistent with the weak decay of the doubly charged baryon 
0. The SELEX result was recently critically discussed in J|. Whether or not the state 
that SELEX reports turns out to be the first observation of doubly charmed baryons, 
studying their properties is important for a full understanding of the strong interaction 
between quarks. 
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Estimations for the masses and spectra of the baryons containing two or more heavy 
quarks have been considered by many authors [|]]. The purpose of this letter is to present 
a consistent treatment of the masses and wave functions of the light, heavy and doubly 
heavy baryons obtained in a simple approximation within the nonperturbative QCD. 
In Ref. || starting from the QCD Lagrangian and assuming the minimal area for the 
asymptotics of the Wilson loop the Hamiltonian of the 3q system in the rest frame has 
been derived. The methodology of the approach has been reviewed recently By 
using this approach and the hypercentral approximation @ we calculate the ground state 
energies and wave functions of the doubly heavy baryons as three quark systems, with 
the three-body confinement force. As a by-product, we also report the masses and wave 
functions for light and heavy baryons. 

From experimental point of view, a detailed discussion of the excited QQ'q states is 
probably premature. Therefore we consider the ground state baryons without radial and 
orbital excitations in which case tensor and spin-orbit forces do not contribute perturba- 
tively. Then only the spin-spin interaction survives in the perturbative approximation. 
The EH has the following form B 



n = VA^ + ^)+Ho + v, (i) 




where H is the non-relativistic kinetic energy operator, V is the sum of the perturbative 
one-gluon exchange potentials V c : 

K (2) 

3 Tin 

and the string potential 

K5tring(n, r 2 , r 3 ) = cr/ min , (3) 

where Z m j n is the sum of the three distances |r»| from the string junction point. In contrast 
to the standard approach of the constituent quark model the dynamical masses are 
no longer free parameters. They are expressed in terms of the running masses mf {Q 2 ) 
defined at the appropriate hadronic scale of Q 2 from the condition of the minimum of the 
baryon mass as function of m,: 



dm ' B ^ 2m, 2 



MF = E T~ + -T + £oK,m 2 ,m 3 ), (4) 



E being eigenvalue of the operator H + V. Technically, this has been done using the 
einbein (auxiliary fields) approach, which is proven to be rather accurate in various calcu- 
lations for relativistic systems. Einbeins are treated as c number variational parameters: 
the eigenvalues of the EH are minimized with respect to einbeins to obtain the physical 
spectrum. Such procedure provides the reasonable accuracy for the meson ground states 



The physical mass Mb of a baryon is |J 

m b = mP + c, c = --T-i (5) 
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where the constant C has the meaning of the quark self energy. The values of r/i are taken 
from 0. They are 1, 0.88, 0.234, and 0.052 for q, s, c, and b quarks, respectively. 

The Effective Hamiltonian is solved using the hyperspherical approach adequate for 
confining potentials. The baryon wave function depends on the three-body Jacobi coor- 
dinates 

IJiij, , x lfJ>ij,k ( rriiVi + mjVj 

Pv = \ — ( Ti ~ r i)' X » = \ — , — ~ Tk 

V /i V A* V mi + rrij 

k cyclic), where ^ and are the appropriate reduced masses 



(6) 



rriimj 
rrii + rrij 



(rrii + m j )m k 
mi + nij + m k 



(7) 



and /i is an arbitrary parameter with the dimension of mass which drops off in the final 
expressions. In terms of the Jacobi coordinates the kinetic energy operator Ho is written 

as 

1 / d 2 d 2 \ 1 f d 2 ^_d_^K 2 {Q)' 
2jl [dp 2 + d>?)~~2^ [M 2 4 



#0 = "^ 



+ RdR 4 



R 2 



where R is the six-dimensional hyper-radius R 2 = p 2 j + \p and K 2 (fl) is angular mo- 
mentum operator whose eigen functions (the hyper spherical harmonics) are K 2 (Q)Y[k] = 
—K(K + 4)yjx], with K being the grand orbital momentum. In terms of Y^] the wave 
function ip(p,\) can be written in a symbolical shorthand as [ID 



ii{ P ,X) = Y,i>K{R)Y [K] (tt). 



(9) 



K 



In the hyper radial approximation which we shall use below K — and ip = ip(R). Since 
R 2 is exchange symmetric the baryon wave function is totally symmetric under exchange. 
Introducing the variable x = ^JJiR and averaging the interaction U = V C + V stT i ng over the 
six-dimensional sphere f2 6 one obtains the Schrodinger equation for u{x) = x°l 2 ip( 



x 



d 2 u(x) 
dx 2 



+ 2 



E + 



bx — 



x 
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u(x) = 0, 



(10) 



with the boundary conditions u(x) 
^( 2/ )~ivr- 1 /V 1/4 exp(-|y 3 / 2 ),y 
state eigenvalue and 

2a, 



; 0(x 5 / 2 
(26) ^ 



cLS X 



-> and the asymptotic u(x) ~ 
oo. In Eq. (|T0|) Eq is the ground 



16 x ^ 
37 ' ^ 



i<3 
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V st rmg(ri,r 2 ,r 3 ) 



dQp 



TP 



ID 



R^fji 

The potential Ktrmg^i, r 2, r 3) has rather complicated structure. In the Y-shape, the 
string meet at 120° in order to insure the minimum energy. This shape moves continuously 
to a two-legs configuration where the legs meet at an angle larger than 120°. Let be 
the angle between the line from quark i to quark j and that from quark j to quark k. If 
(fijk are all smaller than 120°, then the equilibrium junction position coincides with the 
so-called Torichelli point of the triangle in which vertices three quarks are situated. In 
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tt/2 



coscp <-1/2 




2tt/3 tt/3 tc 



Figure 1: The four regions in the (9, x) plane corresponding to (pijk > 120° and < 120° 
for the case of equal quark masses 

this case, in terms of the variables x, 9 = arctan(pi 2 /Ai 2 ), and cos% = pn • 1 p\ 2 \v2 
(0 ^ 9 S tt/2, <; x ^ tt) one obtains 

4 in = X 2 COS 2 ^ (12) 

/ (m? - mi?) tan 2 ^ /m 2 — mi /- \ tan 6* 1 \ 

x / 2 2V+ ^cosx + v^sinx + , 

\mim 2 {mi - m 2 ) \m 2 + m x J m Pi2,s J 

where m 2 = m 1 m 2 m 3 /(m 1 + m 2 + m 3 ). For the case m\ = m 2 = m 3 this expression 
coincides with that derived in Ref. []TTJ . If ^-j. > 120°, the lowest energy configuration 
has the junction at the position of quark j and Z min = + rj k , where 



xsm9 x cos 9 m 2 „ 2m . . 

r 12 = — ■ . ri3 = — W — =-tan 6H tanfcosx + 1, (13) 



and 



xcos6> /m 2 „ 2m „ 

= — \ — ^ tan z 9 tanflcosx + l- (14) 

y/JT^ y m A 2 m 2 

The boundaries corresponding to the various regions in the {9, x) plane are: #1(2) (x) = 
arctan(m 1(2) (=Fcosx-sinx/v / 3)M), 9 3 (x) = arctan(m 2 (/(x) + a// 2 ^) + 4k) /2m), /(%) = 
(1 — «) cosx + (1 + k) sinx/v^3, k = m 1 /m2- These boundaries are shown in Fig. 1 for 
the case of equal quark masses. 
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Note that the frequently used approximation flTjJ, |J is to choose the string junction 
point as coinciding with the center-of-mass coordinate. In this case 

v - — ■ - 32 

•^string ® / ^ m y [^fJ>ij,k 



(7 



15tt 



(15) 

This approximation that greatly simplifies the calculations increases the value of 6 in Eq. 
( [TOP by ~ 5% as illustrated by the results of Table 1. 

We first solve Eq. (|j) for the dynamical quark masses m>i retaining only the string 
potential in the effective Hamiltonian (|1|). This procedure is consistent with Ref. ||, but 
different from that of 
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Then we add the perturbative Coulomb potential and solve 
Eq. (|T0|) to obtain the ground state eigenvalues Eq. The baryon masses Mb are then 
obtained from solving Eq. (|5|). 

Table 1. Illustration of the accuracy of the approximation fll5|) . Shown are the values 
by/m q /a given by Eqs. (|TTD and (15) where m q is the light quark mass. 



barion 


Eq. (11) 


Eq. © 


qqq 


1.583 


1.663 


qqs 


1.556 


1.636 


qss 


1.530 


1.608 


qsc 


1.339 


1.417 


qqb 


1.293 


1.384 


qcb 


1.038 


1.101 


qbb 


0.925 


0.975 



Table 2. For various 3q systems in column (1) we display the dynamical quark masses 
given by Eq. ([|), the ground state eigenvalue Eq in Eq. (|I0"D , the baryon masses including 
the self energy correction Eq. @ (all in units of GeV) and the correlation function 7, 
Eq. (in units GeV 3 / 2 ). 



baryon 


nil 


m 2 


m 3 


Eq 


M B 


7 


qqq 


0.362 


0.362 


0.362 


1.392 


1.144 


0.1389 


qqs 


0.367 


0.367 


0.407 


1.362 


1.242 


0.1369 


qss 


0.371 


0.411 


0.411 


1.335 


1.336 


0.1351 


sss 


0.415 


0.415 


0.415 


1.307 


1.426 


0.1333 


qqc 


0.406 


0.406 


1.470 


1.142 


2.464 


0.1241 


qsc 


0.409 


0.448 


1.471 


1.116 


2.542 


0.1228 


ssc 


0.452 


0.452 


1.473 


1.090 


2.621 


0.1214 


qqb 


0.425 


0.425 


4.825 


1.054 


5.823 


0.1201 


qsb 


0.429 


0.469 


4.826 


1.026 


5.903 


0.1188 


ssb 


0.471 


0.471 


4.826 


1.000 


5.975 


0.1177 


qcc 


0.444 


1.494 


1.494 


0.876 


3.659 


0.1143 


sec 


0.485 


1.496 


1.496 


0.851 


3.726 


0.1134 


qcb 


0.465 


1.512 


4.836 


0.753 


6.969 


0.1136 


scb 


0.505 


1.514 


4.837 


0.729 


7.032 


0.1128 


qbb 


0.488 


4.847 


4.847 


0.567 


10.214 


0.1207 


sbb 


0.526 


4.851 


4.851 


0.544 


10.273 


0.1202 



5 



We use the same parameters as in Ref. a = 0.15 GeV (this value has been 



confirmed in a recent lattice study []I4|), a s = 0.39, mf = 0.009 GeV, mf^ = 0.17 GeV, 
= 1.4 GeV, and m£ = 4.8 GeV. In Table 2 for various three-quark states, we give 
the quark masses m ; , the ground state eigenvalues E , and the baryon masses M B . For 
completeness, in the last column we report the values of the integral 

7 = J ^-dx (16) 

in terms of which the quantities Rijk = (A^ 2 / n 2 ) ■ 7 are expressed determining the 
probability to find a quark i at the location of the quark j in a baryon ijk. These 
quantities are of special importance for the lifetime calculations of heavy hadrons. 

Note that there is no good theoretical reason why quark masses need to be the same 
in different baryons. Inspection of Table 2 shows that the masses of the light quarks (u, 
d or s) are increased by ~ 100 MeV when going from the light to heavy baryons. The 
dynamical masses of light quarks m q ~ yfa ~ 400 — 500 MeV qualitatively agree with the 
results of Ref. flH| obtained from the analysis of the heavy-light ground state mesons. 

While studying Table 2 is sufficient to have an appreciation of the accuracy of our pre- 
dictions, few comments should be added. We expect an accuracy of the baryon predictions 
to be ~ 5 — 10% that is partly due to the approximations employed in the derivations 
of the Effective Hamiltonian itself and partly due to the error associated with the 
variational nature of hyperspherical approximation. From this point of view the overall 
agreement with data is quite satisfactory. E.g. we get AN + A) thcory = 1144 MeV 
vs |(iV + A) cxp = 1085 MeV ( a 5% increase in a s would correctly give the N — A 
center of gravity), ±(A + £ + 2S*) = 1242 MeV vs experimental value of 1267 MeV. We 

1 II 

also find S t heor y = 1336 MeV (without hyperfine splitting) vs H C x P = 1315 MeV and 
2 C theory = 2542 Mev vs H c exp = 2584 MeV. On the other hand, our study shows some 
difficulties in reproducing e.g. the f2-hyperon mass. 

In Table 3 we compare the spin-averaged masses (computed without the spin-spin 



term) of the lowest double heavy baryons to the predictions of other models |T5|, |16| 



|17| as well as variational calculations of Ref. |12| for which the center of gravity of non- 
strange baryons and hyperons is essential a free parameter. Most of recent predictions 
were obtained in a light quark- heavy diquark model [15], [fL6 , in which case the spin- 
averaged values are M = |(Mi/ 2 + M 3 / 2 ). Note that the wave function calculated in the 
hyperspherical approximation shows the marginal diquark clustering in the doubly heavy 
baryons. This is principally kinematic effect related to the fact that in this approximation 
the difference between the various mean values in a baryon is due to the factor a/1 / Jkj 
which varies between ^2/rrii for m 8 = rrij and ^/l/m; for -C rnj. In general, in 
spite of the completely different physical picture, we find a reasonable agreement within 
100 MeV between different predictions for the ground state masses of the doubly heavy 
baryons. Our prediction for M ccu is 3.66 GeV with the perturbative hyperfine splitting 
E* cu — E ccu ~ 40 MeV. Note that the mass of 5+ is rather sensitive to the value of the 

,(°) 



running c-quark mass m c |18 



In conclusion, we have have shown that baryon spectroscopy can be unified in a single 
framework of the Effective Hamiltonian which is consisted with QCD. This picture uses 
the stringlike picture of confinement and perturbative one-gluon exchange potential. The 
main advantage of this work is demonstration of the fact that it is possible to describe all 
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the baryons in terms of the only two parameters inherent to QCD, namely a and a s . 

Table 3. Comparison of our predictions for ground state masses (in units of GeV) of 
doubly heavy baryons with other predictions. 



Baryon 


This 
Work 


Ref. m 


Ref. 15] 


Ref. p| 


Ref. 17] 




3.66 


3.69 


3.57 


3.69 


3.70 


^cc 


3.73 


3.86 


3.66 


3.84 


3.80 


2 C 6 


6.97 


6.96 


6.87 


6.96 


6.99 


Qcb 


7.03 


7.13 


6.96 


7.15 


7.07 


2 bb 


10.21 


10.16 


10.12 


10.23 


10.24 


&bb 


10.27 


10.34 


10.19 


10.38 


10.34 
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